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Abstract 

Background: Secondary edentulism (toothlessness) has evolved on multiple occasions in amniotes including 
several mammalian lineages (pangolins, anteaters, baleen whales), birds, and turtles. All edentulous amniote clades 
have evolved from ancestors with enamel-capped teeth. Previous studies have documented the molecular decay of 
tooth-specific genes in edentulous mammals, all of which lost their teeth in the Cenozoic, and birds, which lost 
their teeth in the Cretaceous. By contrast with mammals and birds, tooth loss in turtles occurred in the Jurassic 
(201.6-145.5 Ma), providing an extended time window for tooth gene degradation in this clade. The release of the 
painted turtle and Chinese softshell turtle genomes provides an opportunity to recover the decayed remains of 
tooth-specific genes in Testudines. 

Results: We queried available genomes of Testudines {Chrysemys picta [painted turtle], Pelodiscus sinensis 
[Chinese softshell turtle]), Aves {Anas platyrhynchos [duck], Gallus gallus [chicken], Meleagris gallopavo [turkey], 
Melopsittacus undulatus [budgerigar], Taeniopygia guttata [zebra finch]), and enamelless mammals {Orycteropus afer 
[aardvark], Choloepus hoffmanni [Hoffmann's two-toed sloth], Dasypus novemcinctus [nine-banded armadillo]) for 
remnants of three enamel matrix protein (EMP) genes with putative enamel-specific functions. Remnants of the 
AMBN and ENAM genes were recovered in Chrysemys and retain their original synteny. Remnants of AM EL were 
recovered in both testudines, although there are no shared frameshifts. We also show that there are inactivated 
copies of AMBN, AMEL and ENAM in representatives of divergent avian lineages including Galloanserae, 
Passeriformes, and Psittaciformes, and that there are shared frameshift mutations in all three genes that predate the 
basal split in Neognathae. Among enamelless mammals, all three EMP genes exhibit inactivating mutations in 
Orycteropus and Choloepus. 

Conclusions: Our results highlight the power of combining fossil and genomic evidence to decipher 
macroevolutionary transitions and characterize the functional range of different loci involved in tooth development. 
The fossil record and phylogenetics combine to predict the occurrence of molecular fossils of tooth-specific genes 
in the genomes of edentulous amniotes, and in every case these molecular fossils have been discovered. The 
widespread occurrence of EMP pseudogenes in turtles, birds, and edentulous/enamelless mammals also provides 
compelling evidence that in amniotes, the only unique, non-redundant function of these genes is in enamel 
formation. 
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Background 

Gnathostomes are unique among vertebrates in utilizing 
jaws and teeth for food acquisition and processing. Diver- 
gent gnathostome lineages exhibit a wide array of dental 
modifications in association with diverse dietary speciali- 
zations. Multiple iterations of secondary tooth loss have 
also evolved in gnathostomes, perhaps most famously in 
birds, turtles, and several mammalian lineages including 
baleen whales, pangolins, and anteaters [1-3]. There are 
also mammalian species with enamelless teeth, i.e., sloths, 
armadillos, aardvarks, pygmy and dwarf sperm whales [2]. 

Tooth development is an intricate process that encom- 
passes a complex series of epithelial-mesenchymal inter- 
actions involving growth factors, transcription factors, 
and signal receptors that affect tooth shape, tooth num- 
ber, and cusp number [4-9]. Several hundred genes are 
associated with tooth development [6], including mem- 
bers of the hedgehog (Hh), fibroblast growth factor (Fgf) 
and bone morphogenic protein (Bmp) families that me- 
diate epithelio-mesenchymal signaling interactions [10]. 
Most or all of the genes that are involved in early tooth 
development are thought to be pleiotropic and have add- 
itional functions outside of tooth formation. However, 
there are also genes with putative tooth-specific functions 
that are expressed later in development by dentin-forming 
odontoblasts and/or enamel-forming ameloblasts [11-17]. 
These genes affect the physical properties of teeth includ- 
ing enamel thickness and structure. The durability of teeth, 
in combination with tooth-specific genes that impact the 
physical structure of dentin and enamel, make teeth a 
model system for studying the coevolution of morpho- 
logical change in the fossil record and molecular change in 
the genome. 



Previous studies have shown that the genes encoding 
three enamel matrix proteins (EMPs), enamelin, amelo- 
genin, and ameloblastin, have become pseudogenized in 
one or more edentulous/enamelless mammals by frameshift 
mutations and/or stop codons [1-3] (Table 1). Further, the 
molecular decay of tooth-specific genes in placental mam- 
mal lineages has been shown to parallel the morphological 
degeneration of enamel in the fossil record [1-3]. The re- 
tention of inactivated EMP genes in mammalian lineages is 
not surprising given that all edentulous and enamelless taxa 
appear to have originated in the Cenozoic [2]. 

By contrast, tooth loss in both birds and turtles occurred 
in the Mesozoic, providing a longer time period for the mo- 
lecular decay of tooth-specific genes in these lineages. In 
birds, the presence of a horny beak and gizzard have pre- 
sumably compensated for edentulism in food acquisition 
and processing [22]. Edentulism in the ancestry of modern 
birds occurred in the Cretaceous, between 125 and 65.5 
Ma, and also evolved independently in several lineages of 
extinct birds that are on the stem to Neornithes (crown 
group birds) [22]. Even though tooth loss in Neornithes oc- 
curred in the Cretaceous, the chicken (Gallus gallus) gen- 
ome retains pseudogenized copies of both AMEL and 
ENAM [18,20,21]. 

Edentulism in Testudines (turtles) occurred even earlier, 
at least as far back as the Late Jurassic [23,24]. The oldest 
turtle is Odontochelys semitestacea from the Triassic of 
China (-220 Ma) [25]. Marginal and palatal teeth are both 
present in this taxon [25]. Slightly younger (-210 Ma) is 
Proganochelys quenstedti from the Upper Triassic [23,26]. 
Marginal teeth are absent in Proganochelys, but palatal 
teeth are present. Phylogenetic analyses suggest that Pro- 
ganochelys, like Odontochelys, is a stem testudine [23]. 



Table 1 Pseudogene remnants of enamel matrix protein (EMP) genes in edentulous and enamelless amniotes 

Tooth gene Taxon 

Turtles Birds Mammals 



1 . ENAM 



2. AMBN 



Frameshift mutations in Chrysemys 
exon 9 (Additional file 1) 



Frameshift mutations in 
Chrysemys picta exons 6 and 1 1 
(Additional file 3) 



3. AMEL 



Framshift mutations in Chrysemys 
picta (exon 1) and Pelodiscus 
sinensis (exons 1, 3, 4) (Additional 
file 7) 



Frameshift mutations in exon 9 of Anas, Gallus, 
Meleagris, Melopsittacus, and Taeniopygia 
(Additional file 1; also see Al-Hashimi et al. [18]) 



Frameshift mutations in Anas platyrhynchos 
(exons 6, 8), Gallus gallus (exons 3, 6), and 
Taeniopygia guttata (exon 6); stop codon in 
exon 8 of Melopsittacus undulatus (Additional 
file 3) 



Frameshift mutations in Anas platyrhynchos 
(exons 1, 2, 3, 4), Gallus gallus (exons 1, 2, 4), 
Meleagris gallopavo (exons 1 , 2), Melopsittacus 
undulatus (exons 1, 2), and Taeniopygia guttata 
(exons 2) (Additional file 7; also see Sire et al. 
[20] and Davit-Beal et al. [21]) 



Frameshift mutations in exon 10 (homologous 
to exon 9 of Alligator) of Tubulidentata 
(aardvark), Pholidota (pangolins), Xenarthra 
(anteaters, sloths, armadillos), Mysticeti (baleen 
whales), and Kogiidae (pygmy and dwarf sperm 
whales) [1,2] 

Frameshift mutations in exon 13 of seven baleen 
whales (Mysticeti) [1,3]; premature stop codon in 
exon 13 of Monodon monoceros [19]; multiple 
inactivating mutations in Orycteropus afer (stop 
codon in exon 5, complete deletion of exon 6, 
frameshift deletion in exon 13) (Additional file 5); 
acceptor splice site mutation in intron 2 of 
Choloepus hoffmanni (AG to AT) 

Frameshift mutations in Choloepus hoffmanni 
(exons 2, 6), deletion of 3' end of intron 6 
(including splice site) and first 14 bp of exon 7 
in Dasypus novemcinctus, and two stop codons 
in exon 6 of Orycteropus afer (Additional file 1 1); 
frameshift mutations in exon 6 of five mysticetes 
[1,3] 
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The oldest testudines with unambiguous crown-group af- 
finities, including tooth loss on the palatines, vomer, and 
pterygoids, are from the Late Jurassic [23] . 

To date, remnants of pseudogenized copies of tooth- 
specific genes in Testudines have not been reported. 
Girondot and Sire [27] were unable to amplify fragments 
of amelogenin in turtles with degenerate PCR primers, 
and it remains unclear if vestiges of tooth-specific genes 
are retained in testudine genomes given the antiquity of 
tooth loss in this clade. On the other hand, phylogenomic 
data suggest that Testudines have slower rates of nuclear 
gene evolution than birds and mammals (fig. 1 in [28]). 
The recent release of two testudine genomes, Chrysemys 
picta (painted turtle) and Pelodiscus sinensis (Chinese 
softshell turtle), provides an opportunity to screen for 
remnants of tooth-specific genes in turtles. These taxa, 
both of which belong to Cryptodira (hidden neck turtles), 
index the basal cladogenic event among crown cryptodires 
and leave only Pleurodira (sideneck turtles) unrepre- 
sented among the three oldest crown-testudine lineages. 
Cryptodira and Pleurodira have traditionally been regarded 
as reciprocally monophyletic based on morphology [29-31] 
and large molecular data sets [32,33]. An exception is 
Barley et al. s [32] coalescence analysis, which recovered 
a basal split between softshells plus the pig-nosed turtle 
(Carettochelys insculpta) and Pleurodira plus other Crypto- 
dira. However, Barley et al.s [32] concatenation analysis 
recovered Cryptodira and Pleurodira. Moreover, the most 
complete testudine phylogeny with a variety of outgroups 
[33] supports the monophyly of both Cryptodira and 
Pleurodira. 

Here, we report the results of querying the genomes 
of Chrysemys and Pelodiscus for inactivated remnants of 
three EMP genes, ENAM, AMEL, and AMBN, all of 
which are hypothesized to have tooth-specific functions. 
We also report the results of querying genome sequences 
of five birds (Gallus gallus [chicken], Meleagris gallopavo 
[turkey], Taeniopygia guttata [zebra finch], Anas platyr- 
hynchos [duck], Melopsittacus undulatus [budgerigar]) and 
three enamelless mammals (Orycteropus afer [aardvark], 
Choloepus hoffmanni [Hoffmanns two-toed sloth], Dasypus 
novemcinctus [nine-banded armadillo]) for remnants of 
these EMP genes. The occurrence of pseudogenized rem- 
nants of EMP genes in highly divergent amniotes, including 
turtles, birds, and edentulous/enamelless mammals, would 
provide robust evidence for the hypothesis that the only es- 
sential, non-redundant function of AMEL, AMBN, and 
ENAM is in enamel formation. Moreover, the search for 
shared frameshift mutations in EMP genes provides an op- 
portunity to date the timing of enamel loss in these clades. 

Methods 

The genomes of Chrysemys picta (painted turtle; PreEn- 
sembl) and Pelodiscus sinensis (Chinese softshell turtle; 



Ensembl 68) were queried with BLASTN using crocodylian 
mRNA sequences for three enamel matrix proteins genes: 
Paleosuchus palpebrosus (Cuviers dwarf caiman) AMEL 
(AF095568), Caiman crocodilus (spectacled caiman) AMBN 
(AY043290), and Crocodylus niloticus (Nile crocodile) 
ENAM (GU344683). Intron-exon boundaries of mRNA 
sequences were determined by blasting crocodylian mRNAs 
against the recently released genome sequence of Alligator 
mississippiensis (American alligator). Crocodylian sequences 
were chosen as probes based on recent evidence that pro- 
vides strong support for a sister group relationship between 
Testudines and Archosauria (crocodylians, birds) [27,34]. 
AMBN has 11 exons in Alligator, all of which include 
protein-coding regions. AMEL includes six exons in 
Alligator, five of which (1-5) contain protein-coding 
regions. Finally, ENAM includes nine exons in Alligator, 
eight of which (2-9) include protein-coding regions. Rem- 
nants of EMP genes that were recovered from Chrysemys 
were subsequently used to query the Pelodiscus genome, 
and remnants of EMP genes from Pelodiscus were used to 
query the Chrysemys genome. Five bird genomes {Anas 
platyrhynchos [duck], Gallus gallus [chicken], Meleagris 
gallopavo [turkey], Melopsittacus undulatus [budgerigar], 
Taeniopygia guttata [zebra finch]) were queried with 
crocodylian mRNAs and virtual pseudogene mRNAs for 
AMEL (EU340348) and ENAM (GU198360) from Gallus, 
as well as with positive hits from the initial round of 
screening for the bird genomes. The NCBI genomes of 
three enamelless mammals {Orycteropus afer [aardvark], 
Choloepus hoffmanni [Hoffmanns two-toed sloth], Dasy- 
pus novemcinctus [nine-banded armadillo]) were queried 
with gene sequences from mammals that retain enamel- 
capped teeth (i.e., Bos taurus, Homo sapiens, Echinops 
telfairi, Elephantulus edwardii, Trichechus manatus). 
Sequences were aligned with Se-Al [35]. 

Selection analyses on branches were performed with 
PAML 4.5 [36] to estimate dN/dS values for an assemblage 
of 20 complete or nearly complete mammalian AMBN 
sequences, all of which were derived from Ensembl 68 
{Ailuropoda melanoleuca, Bos taurus, Callithrix jacchus, 
Equus caballus, Oryctolagus cuniculus) or NCBI {Canis 
lupus familiaris [XM_539304], Choloepus hoffmanni [ABV- 
D01066982, TI 1336370099, ABVD01066984, TI 1338793- 
682, TI 1368206053, ABVD01066985, ABVD01066986], 
Dasypus novemcinctus [AAGV03237582], Echinops telfairi 
[AAIY02097297, AAIY02097298], Elephantulus edwardii 
[AMGZ01205754, AMGZ01205755], Homo sapiens [NM_0- 
16519], Loxodonta africana [AAGU03054364, AAGU0305- 
4363], Orycteropus afer [ALYB0 1124783], Pan paniscus 
[XM_001 160377], Pan troglodytes [XM_001 160377], 
Pongo abelii [XM_0028 14835], Procavia capensis [ABRQ- 
01419216, ABRQ01419217, ABRQ01419218, ABRQ014- 
19219, ABRQ01419220], Sus scrofa [NM_214037]). 
Sequences were aligned with Se-Al [35]. Branch analyses 
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were performed with four branch categories: Orycteropus, 
Choloepus, Dasypus, and all other branches. Orycteropus, 
Choloepus and Dasypus were recognized as separate 
branches in the PAML runs based on the hypothesized in- 
dependent loss of enamel in these lineages [2]. PAML ana- 
lyses were conducted with codon frequency models 1, 2, 
and 3. We used a composite species tree based on Meredith 
et al [37] and Springer et al [38]. Approximate divergence 
dates and fossil ages in Figure 1 are taken from the follow- 
ing sources: Amniota and Testudines to Archosauria [39]; 
Pelodiscus to Chrysemys [40]; Odontochelys [25]; Neog- 
nathae [41]; Galloanserae and Gallus to Meleagris [42]; 
Taeniopygia to Melopsittacus [43]; Ichthyornis [44]; Xenar- 
thra [45]; Xenarthra to Orycteropus [37]; Eomaia [46]. 

Results 

ENAM Gene 

Blast searches with Crocodylus niloticus (Nile crocodile) 
ENAM (GU344683) identified almost the entirety of 



exon 9 in Chrysemys picta (Additional file 1). Remnants 
of this exon in C. picta were identified on scaffold 
JH584398, which has not yet been mapped onto a 
chromosome. The ENAM sequence for exon 9 in Chrys- 
emys includes numerous frameshift indels (Figure 1, 
Table 2, Additional files 1, 2), one of which is a SINE in- 
sertion that shares homology with members of the Cry 
family of SINEs in Cryptodira [47] (Additional files 1, 2). 
An additional SINE insertion occurs 3' to the stop 
codon. Remnants of exon 9 in Chrysemys comprise 
-3840 bp after excluding the SINE insertions. Chrysemys 
exon 9 retains a "TAA" stop codon in the same location 
as Crocodylus. By contrast with the results of blast 
searches against the Chrysemys genome, blast searches 
with Crocodylus ENAM resulted in no significant hits to 
the Pelodiscus sinensis genome. Likewise, blast searches 
with the recovered ENAM pseudogene segment from 
Chrysemys did not find significant matches in the Pelo- 
discus genome. 




Meleagris 



Figure 1 Timetree of edentulous and enamelless amniote taxa for which genomic sequences are available (Ensembl, PreEnsembl, 
NCBI) with a mapping of inferred inactivating mutations (frameshifts, premature stop codons, splice site mutations) in three EMP 
genes to particular branches in the tree. Branches with inactivating mutations are marked by vertical bars (pink = ENAM [E], blue = AMBN [B], 
green = AMEL [A]). Extinct stem taxa with teeth {Odontochelys, Ichthyornis, Eomoio) are also shown. Mutations were mapped to branches by Fitch 
parsimony with delayed transformation optimization. Teeth with enamel, teeth without enamel, and edentulism are denoted with black, gray, 
and white-filled circles, respectively, at terminal and internal nodes. In combination with fossil evidence, dN/dS ratios and frameshift mutations in 
ENAM suggest that enamel was lost independently in the common ancestor of Pilosa (sloths and anteaters) and in multiple armadillo lineages, 
including Dasypus [2]. Remnants of AMBN and ENAM were not recovered from the Pelodiscus genome. References for divergence dates and fossil 
ages are provided in Methods. 
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Table 2 Examples of inactivating mutations in enamel matrix protein genes among edentulous and enamelless 
amniotes with genome sequences 



Taxon 


Inactivating mutations 








ENAM 


AMBN 


AMEL 


Cryptodira 


None 


None 


Possible shared frameshift mutations in exon 2 (see AF9) 


Chrysemys 


AF1 Insertions: 65, 271,411,485, 
994-1 138 (Cry), 2526-2529, 2559- 
2571,3268 

AF1 Deletions: 588, 1264-1276, 
3325-3343 


AF3 Insertions: 956-957, 972, 
1051-1052 

AF3 Deletions: 497-498, 605-61 1, 
794-804, 984-985 


AF7 Deletion: 106 


Pelodiscus 


Gene not recovered 


Gene not recovered 


AF7 Deletions: 122-125, 208, 271, 276 


Neognathae 


AF1 Insertions: 1985-1989,2143- 
2144, 

AF1 Deletions: 2634-2637, 3353- 
3758 


AF3 Deletion: 515-524 


AF7 Deletions: 134 (with subsequent deletion of one additional 
bp in Taeniopygia), 242, 709-710 


Galloanserae 


AF1 Deletion: 1189-1192 


None 


AF7 Deletions: 173-174, 282-289 


Anas 


AF1 Insertions: 552-562, 2031- 
2040 

AF1 Deletions: 2626-2627 


AF3 Deletions: 216-217, 226, 
455, 727-728 


AF7 Insertions: 152, 181 
AF7 Deletion: 88-94 


Gallus + 
Meleogris 


AF1 Insertion: 2654 
AF1 Deletion: 3263-3267 


None 


AF7 Insertion: 484 


Gallus 


AF1 Insertions: 512, 930-931 
AF1 Deletions: 126-127, 585-595 


AF3 Insertions: 225, 230 


AF7 Deletion: 275-276 


Meleogris 


AF1 Insertion: 2438-2441 


Gene not recovered 


AF7 Deletion: 359 


Taeniopygia 


AF1 Deletion: 2958 


None 


None 


+ 

Melopsittocus 








Taeniopygia 


AF1 Deletion: 3318 


AF3 Deletions: 578, 612-613 


AF7 Deletions: 378, 407-417 


Melopsittacus 


AF1 Insertions: 231, 2239 

AF1 Deletions: 432-433, 452-453, 
504-505 


AF3 Stop codon: 708-710 


AF7 Insertion: 105 

AF7 Deletions: 91-94, 387 


Dasypus 


See Meredith et al. [2] 


None 


AF1 1 Deletion: 3' end of intron 6, including the acceptor splice 
site, and the first 14 bp of exon 7 (665-678), which includes the 
last amino acid and the stop codon 


Choloepus 


See Meredith et al. [2] 


Splice site mutation: AG to AT 
at intron 2 acceptor site (gb|ti| 
1338556193) 


AF1 1 Deletions: 45-55, 600 

AF1 1 Initiation codon mutation: ATG to ACG (13-15) 






AF1 1 Stop codon: 356-358 


Orycteropus 


See Meredith et al. [2] 


AF5 Deletions: 415-651 (all of 
exon 6), 1241-1253 

AF5 Stop codon: 388-390 


AF1 1 Stop codons: 239-241, 299-301 



Numbers correspond to sequence alignment positions in Additional files. Additional file 1 = AF1; Additional file 3 = AF3; Additional file 5 = AF5; 
Additional file 7 = AF7; Additional file 9 = AF9; Additional file 11 = AF1 1 . Cry = Cry SINE [47]. 



Blast searches with the virtual pseudogene mRNA of 
Gallus ENAM exon 9 recovered homologous fragments 
of this exon in Meleagris gallopavo (turkey), Taeniopygia 
guttata (zebra finch), Anas platyrhynchos (duck), and 
Melopsittacus undulatus (budgerigar). The 3' region of 
M. gallopavo ENAM was recovered on chromosome 1, 
whereas the 5' region was recovered on the Z chromo- 
some of this species. All of the exon 9 sequences were 
characterized by frameshift mutations (Figure 1), including 



a 406 bp deletion near the 3' end of the coding sequence 
of exon 9 that is shared by all five birds (Additional file 1). 

ENAM sequences for Orycteropus afer, Choloepus 
hoffmanni and Dasypus novemcinctus were reported by 
Meredith et al. [2] and are not duplicated here. Orycteropus 
ENAM includes three frameshift mutations [2]. Among 
xenarthrans, Choloepus ENAM includes numerous frame- 
shift mutations whereas Dasypus ENAM includes a single 
frameshift mutation that is located near the carboxy- 
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terminal end of the coding sequence in exon 10 (= exon 9 
of some non-mammalian vertebrates) [2]. 



AMBN Gene 

Blast searches with Caiman crocodilus (spectacled caiman) 
AMBN (AY043290) recovered portions of exons 1, 2, 3, 6, 
8, 9, 10, and 11, in this order, on Chrysemys scaffold 
JH584398 (Additional files 3 and 4). The AMBN exons are 
located immediately upstream of ENAM exon 9 on scaf- 
fold JH584398. Frameshift mutations occur in exons 6, 10, 
and 11 (Figure 1, Table 2, Additional files 3, 4). AMBN 
sequences were not found in Pelodiscus. 

Among birds, remnants of AMBN were discovered in 
Anas (scaffold 247: exons 2, 3, 5, 6, 7, 8), Gallus (Z 
chromosome: exons 2, 3, 5, 6), Taeniopygia (Z chromo- 
some; exons 6, 7, 8), and Melopsittacus (scaffold JH556633: 
exons 7, 8). Frameshift mutations occur in Anas (exons 2, 
6, 8), Gallus (exon 3) and Taeniopygia (exon 6) (Additional 
file 3). A ten-bp frameshift deletion is shared by Anas, 
Gallus, and Taeniopygia, which together index the com- 
mon ancestry of Neognathae, although this region was not 
recovered in Meleagris and Melopsittacus. There were no 
frameshifts in Melopsittacus exons 7 and 8, but a stop 
codon occurs in the latter exon. 

Among mammals, we recovered exons 1-5, 7, and 10-13 
in Orycteropus afer (Additional file 5). Exon 5 contains a 
stop codon and there is a 13-bp frameshift deletion in exon 
13 (Table 2). Exon 6 has been deleted along with portions 
of introns 5 and 6 (Table 2). Exons 8 and 9 are 39-bp dupli- 
cations of exon 7 that occur in primates [48] and their ab- 
sence in Orycteropus is not unexpected. With the exception 
of Elephantulus edwardii, which has two duplications of 
exon 7 as in Homo, other afrotherians (Loxodonta africana, 
Trichechus manatus, Procavia capensis, Chrysochloris asia- 
tica, Echinops telfairi) are similar to Orycteropus in posses- 
sing only exon 7. We also recovered sequences for exons 
1-7 and 10-13 in Dasypus AMBN (Additional file 5). The 
recovered protein-coding sequence is intact, although the 
presumed start codon is ten codons downstream of the 
start codon in the human sequence and occurs in exon 2 
rather than exon 1, as is also the case for Loxodonta afri- 
cana (African elephant) (Additional file 6). Complete cod- 
ing sequences for nine exons (1-5, 7, 10-12) and partial 
sequences for exons 6 and 13 were recovered from the Cho- 
loepus genome (Additional file 5). The putative start codon 
in Choloepus is located in exon 2 as in Dasypus. There are 
no frameshift mutations or stop codons in the available 
Choloepus sequence, but intron 2 exhibits an inactivating 
mutation at the acceptor splice site (AG to AT; see TI# 
1338556193 in Trace Archives) (Table 2). By contrast, all of 
the donor and accepter splice sites exhibit canonical GT 
and AG motifs, respectively, in 16 mammals with enamel- 
capped teeth (taxon names provided in Methods). 



Selection analyses based on an alignment of 20 mamma- 
lian AMBN sequences (Additional file 6) with codon fre- 
quency (CF) models 1, 2, and 3 all suggest elevated dN/dS 
ratios on the Orycteropus (CF1 = 0.71, CF2 = 0.77, CF3 = 
0.89), Choloepus (CF1 = 0.86, CF2 = 0.90, CF3 = 1.04) and 
Dasypus (CF1 = 0.64, CF2 = 0.71, CF3 = 0.74) branches 
relative to the background rate (CF1 = 0.48, CF2 = 0.51, 
CF3 = 0.56) in mammalian taxa with enamel-covered 
teeth. 

AMEL Gene 

Blast searches with Paleosuchus palpebrosus (Cuviers 
dwarf caiman) AMEL (AF095568) recovered most of 
exons 1, 3, 4, and 5 in Pelodiscus, with ancestral synteny 
preserved, on scaffold JH208023 (Additional files 7 and 8). 
Exons 1 and 5 were identified on scaffold JH584884 in 
Chrysemys (Additional files 7 and 8). Chrysemys and Pelo- 
discus both retain the same start codon as in crocodylians. 
Frameshift mutations are present in sequences of both tes- 
tudine species (exon 1 in Chrysemys, exons 1, 3, and 4 in 
Pelodiscus), although there are no shared frameshifts in 
these two exons (Figure 1, Table 2, Additional files 7 and 
8). In addition, Pelodiscus and Chrysemys share exon- 
flanking sequences that are 3' to exon 1 and 5' to exon 5. 
The sequence that flanks exon 1 is > 1000 bp and includes 
regions that are homologous to intron 1 and possibly exon 
2 of Alligator mississippiensis (Additional file 6). Intron 1 
retains the canonical "GT" splice donor site in both testu- 
dines (Additional file 9). The homology of exon 2 is less 
certain, but nevertheless of potential importance as there 
are multiple shared frameshift mutations in Chrysemys 
and Pelodiscus (Additional file 9). The sequence that 
flanks exon 5 is > 600 bp long and is homologous to a por- 
tion of intron 4 of A. mississippiensis with the canonical 
"AG" splice acceptor site in Chrysemys ("AA" in Pelodis- 
cus) (Additional file 10). 

In birds, blast searches with Gallus AMEL pseudogene 
mRNA (EU340348) recovered homologous DNA sequences 
in Anas (exons 1-4 with original synteny on scaffold 415), 
Meleagris (exons 1-4 with original synteny on chromosome 
1), Melopsittacus (exons 1-4 with original synteny on 
AGAI01063440), and Taeniopygia (exons 2-4 with original 
synteny on chromosome 1) (Additional file 7). All of the 
birds except for Anas share a common start codon with 
crocodylians. Frameshift mutations occur in all avian taxa, 
including a one-bp deletion in exon 4 that is shared by all 
five birds (Figure 1). 

Among mammals, there are two stop codons in exon 
6 of Orycteropus AMEL (Figure 1, Table 2, Additional 
file 11). Choloepus AMEL shows replacement of the 
methionine initiation codon in exon 2 by a threonine 
codon, an 11 bp frameshift deletion in exon 2, and a 
premature stop codon followed by a frameshift deletion 
in exon 6 (Figure 1, Table 2, Additional file 11). All of 
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the inactivating mutations in Choloepus are corroborated 
by chromatograms in NCBIs Trace Archives. The se- 
quence for Dasypus is intact through exon 6 (excepting 
exon 4, see below), but a deletion incorporates the 3' end 
of intron 6, including the acceptor splice site, and the first 
14 bp of exon 7, including the stop codon. Orycteropus, 
Choloepus, and Dasypus also lack a functional copy of 
exon 4, either because of stop codons (Orycteropus, Cho- 
loepus) or a donor splice site mutation (GT to AT) in in- 
tron 4 that is adjacent to this exon. However, a functional 
copy of exon 4 is variably present in Mammalia [48] and is 
missing from six other afrotherians (Chrysochloris, Echi- 
nops, Elephantulus, Loxodonta, Procavia, Trichechus) with 
genome sequences because of stop codons, frameshift 
mutations, and/or splice site mutations (data not shown). 

Discussion 

EMP Genes in Turtles 

Crown Testudines comprise 319 extant species [49], all 
of which are toothless and instead have a keratinized 
beak that in combination with strong jaw muscles allows 
testudines to tear food and capture prey [21]. The oldest 
testudine fossil is the Late Triassic Odontochelys semites- 
tacea from China [25]. Odontochelys lacked a beak and 
instead retained both marginal and palatal teeth [25]. 
Proganochelys quenstedti is also known from the Late 
Triassic, and is more derived than Odontochelys. Proga- 
nochelys lacked marginal teeth, but palatal teeth were 
present. The morphology of the mandibles suggests that 
Proganochelys possessed a keratinized beak [26] . It remains 
unclear if palatal teeth, which occur on the vomer, palatine, 
and pterygoid bones, were lost in the common ancestor of 
all living turtles or independently in cryptodires and pleur- 
odires owing to alternate phylogenetic hypotheses wherein 
taxa with palatal teeth such as Kayentachelys apix are posi- 
tioned as stem testudines or as stem cryptodires [23,24,29- 
31,50,51]. In either case palatal teeth were lost no later 
than the Late Jurassic in the ancestry of crown testudines. 

The retention of pseudogenized copies of three EMP 
genes (AMEL, AMBN, ENAM) in Chrysemys is perhaps 
surprising given the antiquity of tooth loss in turtles and 
the extended time window for large-scale deletions and/or 
rearrangements to erase or scramble the genomic instruc- 
tions for enamel production in this edentulous taxon. 
However, remnants of AMBN and ENAM exon 9 in 
Chrysemys occur in juxtaposition to each other on the 
same contig. Kawasaki's [52] reconstructions of SCPP gene 
order in both stem tetrapods and stem amniotes suggest 
that AMBN is immediately upstream of ENAM. Thus, it 
appears that remnants of this original gene order are 
present in the painted turtle genome. The protein-coding 
regions of AMBN and ENAM in Chrysemys have been bat- 
tered by inactivating mutations, including SINE insertions 
in exon 9 of ENAM, but nevertheless retain unambiguous 



signatures of their heritage. The only EMP gene that was 
recovered in Pelodiscus is AMEL. It remains unclear 
whether AMBN and ENAM were completely deleted from 
the Pelodiscus genome during evolutionary history, are dif- 
ficult to recognize because of numerous mutations, or are 
missing from the current assembly of the Pelodiscus gen- 
ome due to incomplete sequencing assembly. Pelodiscus 
AMEL is relatively complete and includes exons 1, 3, 4, 
and 5 in their original syntenic order. By contrast, we only 
discovered exons 1 and 5 of Chrysemys AMEL. The two 
exons represent the only overlapping EMP fragments that 
were discovered in both Pelodiscus and Chrysemys. Frame- 
shift mutations in AMEL occur in both species, but the 
only shared frameshift mutations occur in putative exon 2, 
which has equivocal homology with the same exon in cro- 
codylians. Nevertheless, if these frameshift mutations are 
genuine then there is molecular evidence for the inactiva- 
tion of enamel production in the common ancestry of 
Cryptodira. 

The genetic and developmental basis of edentulism in 
turtles involves both tooth-specific and pleiotropic genes. 
Results presented here demonstrate that three tooth- 
specific EMP genes were pseudogenized in the ancestry of 
modern cryptodires. Previously, Tokita et al. [53] exam- 
ined patterns of gene expression in the developing lower 
jaw of Pelodiscus sinensis and reported impairment of Shh 
signaling in the oral epithelium along with early-stage ar- 
rest of odontoblast development by abrogation of Msx2 
expression in dental mesenchyme. It remains unclear if 
changes in gene expression or EMP pseudogenization oc- 
curred first in the evolutionary history of turtles. 

EMP Genes in Birds 

AMEL and ENAM pseudogenes have previously been 
reported in Gallus [20], which belongs to Galliformes. 
We provide extended evidence for inactivated EMP 
genes in birds and show that there are remnants of 
AMEL and ENAM in four other birds with complete gen- 
ome sequences: Meleagris (Galliformes), Anas (Anseri- 
formes), Melopsittacus (Psittaciformes), and Taeniopygia 
(Passeriformes) (Table 1). In addition, we recovered mo- 
lecular evidence for an inactivated copy of AMBN in 
Anas, Gallus, Taeniopygia, and Melopsittacus. It is note- 
worthy that there are putative frameshift mutations in 
AMBN, AMEL, and ENAM that are shared by represen- 
tatives of Galliformes (Gallus, Meleagris), Anseriformes 
(Anas), and Neoaves (Taeniopygia, Melopsittacus), which 
together index the deepest split in Neognathae [54]. It 
will be important to determine if these inactivating muta- 
tions are shared with Palaeognathae (tinamous and 
ratites), in which case molecular evidence would be con- 
sistent with the loss of enamel in the common ancestor 
of Neornithes rather than independently in Neognathae 
and Palaeognathae. 
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EMP Genes in Mammals 

It has previously been shown that ENAM is a pseudogene 
in a wide array of edentulous and enamelless mammals 
[1,2] (Table 1). Further, AMEL and AMBN sequences have 
been reported for baleen whales and in both cases are 
inactivated in multiple species [1,3]. However, the func- 
tional versus pseudogene status of AMEL and AMBN 
remains to be investigated in most edentulous and enam- 
elless taxa. 

Here, we provide evidence of inactivating mutations in 
AMEL in Orycteropus afer (aardvark) and Choloepus 
hoffmanni (Hoffmanns two-toed sloth), both of which 
have enamelless teeth comprised of dentin. Multiple 
frameshift mutations in Choloepus AMEL are confirmed 
by Trace Archives chromatograms (Additional file 8). 
The two stop codons in Orycteropus are based on Illu- 
mini Hi-Seq sequencing technology with 44X genome 
coverage. 

Previous studies suggest that AMBN may play a role 
in dentin formation/regeneration [55,56]. However, this 
hypothesis is contradicted by our finding that there are 
inactivating mutations in Orycteropus and Choloepus 
AMBN, Moreover, dN/dS ratios indicate relaxed purify- 
ing selection in these taxa. Similarly, McGowen [19] 
reported a premature stop codon in exon 13 (JF504758) 
of Monodon monoceros (narwhal). Narwhals have ves- 
tigial teeth, and in the case of males a single enlarged 
tusk. Both the vestigial teeth and tusk are composed of 
dentin and cementum, but no enamel [57-59]. The oc- 
currence of a premature stop codon in Monodon AMBN, 
in conjunction with inactivating mutations in two other 
lineages with enamelless teeth {Orycteropus, Choloepus), 
is consistent with the hypothesis that the only essential, 
non-redundant role of AMBN is in enamel formation. 
Along these lines, ameloblastin knockout mice do not 
form an enamel layer, and lack dentin defects [13,15]. 

By contrast with Choloepus, where enamel loss occurred 
more than 50 million years ago, enamel degeneration in 
Dasypus occurred more recently [2]. Dasypus retains a 
vestigial enamel or enamel-like substance that covers the 
dentin of some teeth and is quickly worn off. Moreover, 
there are no frameshifts in ENAM that are shared by 
Dasypus and other armadillo genera [2], The single frame- 
shift in Dasypus ENAM precludes a functional, full-length 
enamelin protein, but the location of the frameshift is 
close to the 3' end of the protein-coding sequence in exon 
10, and shorter protein products are possible [2]. Our 
finding that AMBN is intact in Dasypus, and that only the 
last amino acid is truncated from the coding region of 
AMEL, is consistent with the presence of vestigial enamel' 
in this species [2]. We also note that a partial AMBN se- 
quence (JF701624) for Kogia breviceps (pygmy sperm 
whale) has an intact protein-coding region even though 
K breviceps has enamelless teeth and frameshift mutations 



in the ENAM gene. It will be important in future studies 
to sequence the remaining exons of Kogia AMBN to de- 
termine if the missing exonic regions contain frameshift 
mutations and/or other inactivating mutations. 

Conclusions 

Edentulism has evolved independently in multiple lineages 
of living amniotes including turtles, birds, echidnas, baleen 
whales, anteaters, and pangolins. There are also mammals 
with enamelless teeth including pygmy sperm whale, nar- 
whal, sloths, armadillos, and aardvarks. In every case these 
edentulous or enamelless forms have descended from 
ancestors with enamel-capped teeth. Thus, amniote di- 
versity provides a natural laboratory for testing hypoth- 
eses of tooth-specific gene function [3]. Moreover, this 
laboratory includes multiple, replicated experiments. 
AMBN, AMEL, and ENAM have all been postulated to 
have tooth-specific or even enamel-specific gene func- 
tions [11-17], although pleiotropic functions have been 
suggested for AMBN [9,48,60,61] and AMEL [62,63]. The 
widespread occurrence of EMP pseudogenes in turtles, 
birds, and several mammalian lineages (Table 1) provides 
compelling evidence that the only unique, non-redundant 
function of these genes in amniotes is in enamel forma- 
tion: functional copies of these genes have not been 
retained by natural selection when enamel production was 
abrogated independently in distantly related lineages. By 
contrast, representative sequences from amniotes with 
enamel-capped teeth retain intact coding sequences for 
AMBN [1,48]; also see Additional file 11, AMEL [3,48,64], 
and ENAM [1,2,18,37,65]. 

The evolution of tooth loss in multiple amniote lineages 
also provides a model system for integrating the fossil rec- 
ord, phylogenetics, and genomics. This system allows for 
reciprocal hypothesis testing and provides a multifaceted, 
synthetic view on macroevolutionary transitions in testu- 
dines, birds, and edentulous/enamelless mammals [1,2]. 
The fossil record and phylogenetics combine to predict 
the occurrence of molecular fossils of tooth-specific genes 
in the genomes of edentulous and enamelless amniotes, 
and molecular fossils have been discovered in all lineages 
that have been investigated. We suggest that these striking 
patterns may prove especially useful for educating the 
public on the convergence of evidence that clearly docu- 
ments the evolutionary process over deep time. 

Additional files 



Additional file 1 : Nexus alignment of ENAM exon 9 sequences for 
two crocodylians, one testudine, and five birds. The stop codon in 
crocodylians and Chrysemys picta occurs at positions 3788-3790. A 406 bp 
deletion that is shared by all five birds occurs near the 3' end of the coding 
region at positions 3353-3758. Additional frameshift mutations are shared 
by four of five birds including representatives of Anseriformes, Galliformes, 
and Passeriformes (2634-2637). Cry SINE elements [47] in Chrysemys occur at 
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positions 994-1 138 and 4782-4994. Several long strings of Ns in Chrysemys 
picta were abridged to shorter strings often Ns. 

Additional file 2: Schematic alignment showing regions of ENAM 
exon 9 that were identified in Chrysemys picta. Green bars below 
sequence numbers show regions of sequence similarity between C. picta 
and Alligator mississippiensis. Red rectangles = frameshift mutations in C. 
picta; light blue rectangles = indels within the coding sequence that are 
in multiples of three base pairs; CRY = Cry SINE insertion; red arrow = 
position of stop codon in A. mississippiensis. 

Additional file 3: Nexus alignment of AMBN exon sequences for 
two crocodylians, one testudine, and one bird. The start and stop 
codons in crocodylians occur at positions 127-129 (exon 2) and 1374- 
1376 (exon 11), respectively. 

Additional file 4: Schematic alignment showing regions of AMBN 
exons 1-11 that were identified in Chrysemys picta. Green bars below 
sequence numbers show regions of sequence similarity between C. picta 
and Alligator mississippiensis. Dark blue rectangles = numbered exons in 
A. mississippiensis; red rectangles = frameshift mutations in C. picta; light 
blue rectangles = indels within the coding sequence that are in multiples 
of three bp; dark blue arrow = position of start codon in A. 
mississippiensis; red arrow = position of stop codon in A. mississippiensis. 

Additional file 5: Nexus alignment of AMBN exon sequences for 
Homo sapiens, Orycteropus afer and six other afrotherians, Dasypus 
novemcinctus, and Choloepus hoffmanni. The start and stop codons in 
Homo occur at positions 121-123 (exon 1) and 1468-1470 (exon 13), 
respectively. The presumed start codon in Dasypus occurs at positions 151- 
153. Dasypus shares the same stop codon as Homo at positions 1468-1470, 
whereas the presumed stop codon in Choloepus is at positions 1477-1479. 

Additional file 6: Nexus alignment of protein-coding AMBN 
sequences for 20 mammals that was used in selection analyses with 
PAML 4.5 [36]. 

Additional file 7: Nexus alignment of AM EL exon sequences in nexus 
format for two crocodylians, two testudines, and five birds. The start 
and stop codons in crocodylians occur at positions 69-71 (exon 1) and 731- 
733 (exon 5), respectively. A putative one-bp frameshift deletion that is 
shared by all five birds species occurs at position 242 in exon 4. 

Additional file 8: Schematic alignment showing regions of AMEL 
exons 1 -5 that were identified in Chrysemys picta and Pelodiscus 
sinensis. Bars below sequence numbers show regions of sequence similarity 
between Alligator mississippiensis and both testudines (green) or A 
mississippiensis and one testudine (chartreuse). Dark blue rectangles = 
numbered exons in A mississippiensis; red rectangles = frameshift mutations 
in C. picta; light blue rectangles = indels within the coding sequence that 
are in multiples of three bp; dark blue arrow = position of start codon in A 
mississippiensis; red arrow = position of stop codon in A mississippiensis. 

Additional file 9: Nexus alignment of Alligator mississippiensis AMEL 
exon 1 (positions 1-114), intron 1 (positions 115-813), and exon 2 
(838-894) with exon 1 and flanking sequences from Chrysemys 
picta and Pelodiscus sinensis. 

Additional file 10: Nexus alignment of AMEL intron 4 (positions 1- 
662) and exon 5 (positions 663-844) sequences for Chrysemys picta 
and Pelodiscus sinensis. 

Additional file 11: Nexus alignment of AMEL sequences (exons 2-7) 
for Homo sapiens, Orycteropus afer and six other afrotherians 
{Echinops telfairi, Chrysochloris asiatica, Elephantulus edwardii, 
Trichechus manatus, Loxodonta africana, Procavia capensis), Dasypus 
novemcinctus, and Choloepus hoffmanni. The start and stop codons in 
Homo occur at positions 13-15 (exon 2) and 668-670 (exon 7), respectively. 
The mutated stop codon and frameshift mutation in exon 2 of Choloepus 
are clear in gnl|ti|1 3252891 35 and gnl|ti| 131 1883086. The stop codon in 
exon 4 of Choloepus is confirmed by traces gnl|ti|1 325336699 and gnl|ti| 
1312095734. The frameshift mutation in exon 6 of Choloepus is evident in 
trace gnl|ti| 13635791 07. The deletion that includes the 3' region of intron 6 
and first 14 bp of exon 7 in Dasypus is evident in gn|ti| 1886397515. 



Competing interests 

The authors declare that they have no competing financial interests. 



Authors' contributions 

JG, RWM, and MSS conceived and designed the study. RWM and MSS 
performed the analyses. MSS wrote the manuscript. JG and RWM provided 
comments on the manuscript. JG and MSS prepared figures. All authors read 
and approved the final manuscript. 

Acknowledgements 

This work was supported by NSF (EF0629860 to MSS and JG; DEB0743724 to 
JG). We thank four anonymous referees for helpful comments on an earlier 
draft of this manuscript. 

Author details 

department of Biology, University of California, Riverside, CA 92521, USA. 
2 Current address: Department of Biology and Molecular Biology, Montclair 
State University, Montclair, NJ 07043, USA. 

Received: 30 October 2012 Accepted: 16 January 2013 
Published: 23 January 2013 

References 

1. Demere TA, McGowen MR, Berta A, Gatesy J: Morphological and molecular 
evidence for a stepwise evolutionary transition from teeth to baleen in 
mysticete whales. Syst Biol 2008, 57:15-37. 

2. Meredith RW, Gatesy J, Murphy WJ, Ryder OA, Springer MS: Molecular 
decay of the tooth gene enamelin {ENAM) mirrors the loss of enamel in 
the fossil record of placental mammals. PLoS Genet 2009, 5:e1 000634. 

3. Meredith RW, Gatesy J, Cheng J, Springer MS: Pseudogenization of the 
tooth gene enamelysin {MMP20) in the common ancestor of extant 
baleen whales. Proc R Soc B 201 1, 278:993-1002. 

4. Scarel-Caminaga RM, Pasetto S, da Silva ER, Peres RCR: Genes and tooth 
development: reviewing the structure and function of some key players. 
BrazJOral Sci 2003, 2:339-347. 

5. Chai Y, Maxson RE Jr: Recent advances in craniofacial morphogenesis. Dev 
Dynamics 2006, 235:2353-2375. 

6. Thesleff I: The genetic basis of tooth development and dental defects. 
Am J Med Gen 2006, 140A:2530-2535. 

7. Kapadia H, Mues G, D'Sousa R: Genes affecting tooth morphogenesis. 
Orthod Craniofacial Res 2007, 10:105-1 13. 

8. Bei M: Molecular genetics of tooth development. Curr Opin Gen Dev 2009, 
19:504-510. 

9. Caton J, Tucker AS: Current knowledge of tooth development: patterning 
and mineralization of the murine dentition. J Anat 2009, 214:502-515. 

10. Fleischmannova J, Matalova E, Tucker AS, Sharpe PT: Mouse models of 
tooth abnormalities. Eur J Oral Sci 2008, 1 16:1-10. 

11. Snead ML, Lau EC, Zeichner-David M, Fincham AG, Woo SLC, Slavkin HQ 
DNA sequence for cloned cDNA for murine amelogenin reveal the 
amino acid sequence for enamel-specific protein. Biochem Biophys Res 
Comm 1985, 129:812-818. 

12. Rajpar MH, Harley K, Laing C, Davies RM, Dixon MJ: Mutation of the gene 
encoding the enamel-specific protein, enamelin, causes autosomal- 
dominant amelogenesis imperfecta. Hum Mol Gen 2001, 10:1673-1677. 

13. Fukumoto S, Kiba T, Hall B, lehara N, Nakamura T, Longenecker G, Krebsbach 
PH, Nanci A, Kulkarni AB, Yamada Y: Ameloblastin is a cell adhesion 
molecule required for maintaining the differentiation state of 
ameloblasts. J Cell Biol 2004, 167:973-983. 

14. Fukumoto S, Arakaki M, Iwamoto T, Yamada A, Miyamoto R, Naruse M, 
Nakamura T: Epithelial cell lines in the field of dental research: review. In 
Interface of Oral Health Science 201 1. Edited by Sasaki K, Suzuki O, Takahashi 
N. New York: Springer; 2012:327-333. 

15. Hu JC-C, Hu Y, Smith CE, McKee MD, Wright JT, Yamakoshi Y, Papagerakis P, 
Hunter GK, Feng JQ, Yamakoshi F, Simmer JP: Enamel defects and 
ameloblast-specific expression in Enam knock-out//c?cZ knock-in mice. 

J Biol Chem 2008, 283:10858-10871. 

16. Hatakeyama J, Fukumoto S, Nakamura T, Haruyama N, Suzuki S, Hatakeyama 
Y, Shum L, Gibson CW, Yamada Y, Kulkarni AB: Synergistic roles of 
amelogenin and ameloblastin. J Dent Res 2009, 88:318-322. 

17. Chan H-C, Estrella NMRP, Milkovich RN, Kim J-W, Simmer JP, Hu JC-C: Target 
gene analyses of 39 amelogenesis imperfecta kindreds. Eur J Oral Sci 
2011, 119(Suppl 1)311-323. 

18. Al-Hashimi N, Lafont A-G, Delgado S, Kawasaki K, Sire J-Y: The enamelin 
genes in lizard, crocodile, and frog and the pseudogene in the chicken 



Meredith et al. BMC Evolutionary Biology 2013, 13:20 
http://www.biomedcentral.eom/1 47 1 -2 1 48/1 3/20 



Page 10 of 10 



provide new insights on enamelin evolution in tetrapods. Mol Biol Evol 
2010, 27:2078-2094. 

19. McGowen MR: Toward the resolution of an explosive radiation - a 
multilocus phylogeny of oceanic dolphins (Delphinidae). Mol Phylogenet 
Evol 2011,60:345-357. 

20. Sire J-Y, Delgado SC, Girondot M: Hen's teeth with enamel cap: from 
dream to impossibility. BMC Evol Biol 2008, 8:246. 

21 . Davit-Beal T, Tucker AS, Sire J-Y: Loss of teeth and enamel in tetrapods: 
fossil record, genetic data and morphological adaptations. J Anat 2009, 
214:477-501. 

22. Louchart A, Viriot L: From snout to beak: the loss of teeth in birds. Trends 
EcolEvol 2011, 26:663-673. 

23. Joyce WG: Phylogenetic relationships of Mesozoic turtles. Bull Peabody 
Mus Nat Hist 2007, 48:3-102. 

24. Anquetin J, Barrett PM, Jones MEH, Moore-Fay S, Evans SE: A new stem turtle 
from the Middle Jurassic of Scotland: new insights into the evolution and 
palaeoecology of basal turtles. Proc Roy Soc B 2009, 276:879-886. 

25. Li C, Wu X-C, Rieppel O, Wang L-T, Zhao L-J: An ancestral turtle from the 
Late Triassic of southwestern China. Nature 2008, 456:497-501. 

26. Gaffney ES: The comparative osteology of the Triassic turtle 
Progonochelys. Bull Amer Mus Nat Hist 1 990, 1 94:1 -263. 

27. Girondot M, Sire J-Y: Evolution of the amelogenin gene in toothed and 
toothless vertebrates. Eur J Oral Sci 1998, 106(Suppl 1):501 -508. 

28. Chiari Y, Cahais V, Galtier N, Delsuc F: Phylogenomic analyses support the 
position of turtles as the sister group of birds and crocodiles 
(Archosauria). BMC Biol 2012, 10:65. 

29. Gaffney ES, Hutchison JH, Jenkins FA Jr, Meeker LJ: Modern turtle origins: 
the oldest cryptodire. Science 1987, 237:289-291. 

30. Gaffney ES, Meylan PA, Wyss A: A computer assisted analysis of the 
relationships of the higher categories of turtles. Cladistics 1991, 7:313-335. 

31 . Gaffney ES, Jenkins FA Jr: The cranial morphology of Kayentachelys, an 
Early Jurassic cryptodire, and the early history of turtles. Acta Zool 2010, 
91:335-368. 

32. Barley AJ, Spinks PQ, Thomson RC, Shaffer HB: Fourteen nuclear genes 
provide phylogenetic resolution for difficult nodes in the turtle tree of 
life. Mol Phylogenet Evol 2010, 55:1 189-1 194. 

33. Guillon J-M, Guery L, Hulin V, Girondot M: A large phylogeny of turtles 
(Testudines) using molecular data. Contrib Zool 2012, 81:147-158. 

34. Crawford NG, Faircloth BC, McCormack JE, Brumfield RT, Winker K, Glenn TC: 
More than 1000 ultraconserved elements provide evidence that turtles 
are the sister group of archosaurs. Biol Lett 2012, 8:783-786. 

35. Rambaut A: Se-AI: Sequence Alignment Editor, 2.0a 7 7.1 996. http://evolve.zoo. 
ox.ac.uk. 

36. Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol 
Evol 2007, 24:1586-1591. 

37. Meredith RW, Janecka JE, Gatesy J, Ryder OA, Fisher CA, Teeling EC, 
Goodbla A, Eizirik E, Simao TLL, StadlerT, Rabosky DL, Honeycutt RL, Flynn 
JJ, Ingram CM, Steiner S, Williams TL, Robinson T, Burk-Herrick A, Westerman 
M, Ayoub NA, Springer MS, Murphy WJ: Impacts of the Cretaceous 
Terrestrial Revolution and KPg extinction on mammalian diversification. 
Science 201 1,334:521-524. 

38. Springer MS, Meredith RW, Gatesy J, Emerling CA, Park J, Rabosky DL, 
Stadler T, Steiner C, Ryder OA, Janecka JE, Fisher CA, Murphy WJ: 
Macroevolutionary dynamics and historical biogeography of primate 
diversification inferred from a supermatrix. PLoS ONE 2012, 7:e45921. 

39. Shedlock AM, Edwards SC: Amniotes (Amniota). In The Timetree of Life. 
Edited by Hedges SB, Kumar S. Oxford: Oxford Univ Press; 2009:375-379. 

40. Shaffer HB: Turtles (Testudines). In The Timetree of Life. Edited by Hedges 
SB, Kumar S. Oxford: Oxford Univ Press; 2009:398-401. 

41. Pereira SL, Baker AJ: A mitogenomic timescale for birds detects variable 
phylogenetic rates of molecular evolution and refutes the standard 
molecular clock. Mol Biol Evol 2006, 23:1731-1740. 

42. Pereira SL, Baker AJ: A molecular timescale for galliform birds accounting 
for uncertainty in time estimates and heterogeneity of rates of DNA 
substitutions across lineages and sites. Mol Phylogenet Evol 2006, 
38:499-509. 

43. van Tuinen M: Advanced birds (Neoaves). In The Timetree of Life. Edited by 
Hedges SB, Kumar S. Oxford: Oxford Univ Press; 2009:419-422. 

44. Clarke JA: Morphology, phylogenetic taxonomy, and systematics of 
Ichthyornis and Apatomis (Avialae: Ornithurae). Am Mus Novitat 2006, 
286:1-179. 



45. Delsuc F, Douzery EJP: Armadillos, anteaters, and sloths (Xenarthra). In The 

Timetree of Life. Edited by Hedges SB, Kumar S. Oxford: Oxford Univ Press; 
2009:475-478. 

46. Ji Q, Luo Z-X, Yuan C-X, Wible JR, Zhang J-P, Georgi JA: The earliest known 
eutherian mammal. Nature 2002, 416:816-822. 

47. Sasaki T, Takahashi K, Nikaido M, Miura S, Yasukawa Y, Okada N: First 
application of the SINE (short interspersed repetitive element) method 
to infer phylogenetic relationships in reptiles: an example from the 
turtle superfamily Testudinoidea. Mol Biol Evol 2004, 21:705-715. 

48. Sire J-Y, Davit-Beal T, Delgado S, Gu X: The origin and evolution of enamel 
mineralization genes. Cells Tissue Organs 2007, 186:25-48. 

49. Turtle Taxonomy Working Group: An annotated list of modern turtle 
terminal taxa with comments on areas of taxonomic instability and 
recent change. Chelonian Res Monogr 2007, 4:1 73-1 99. 

50. Sterli J, Joyce WG: The cranial anatomy of the lower Jurassic turtle 
Kayentachelys apix. Acta Palaentol Polon 2007, 52:675-694. 

51 . Sterli J: A new, nearly complete stem turtle from the Jurassic of South 
America with implications for turtle evolution. Biol Lett 2008, 4:286-289. 

52. Kawasaki K: The SCPP gene repertoire in bony vertebrates and graded 
differences in mineralized tissues. Dev Gene Evol 2009, 219:147-157. 

53. Tokita M, Chaeychrmsri W, Siruntawineti J: Developmental basis of 
toothlessness in turtles: insight into convergent evolution of vertebrate 
morphology. Evolution 2012, doi:10.1 1 1 1/j.1 558-5646.201 2.01 752.x. 

54. Hackett SJ, Kimball RT, Reddy S, Bowie RCK, Braun EL, Braun MJ, Chojnowski 
JL, Cox WA, Han K-L, Harshman J, Huddleston CJ, Marks BD, Miglia KJ, 
Moore WS, Sheldon FH, Steadman DW, Witt CC, Yuri T: A phylogenomic 
study of birds reveals their evolutionary history. Science 2008, 
320:1763-1768. 

55. Nakamura Y, Slaby I, Spahr A, Pezeshki G, Matsumoto K, Lyngstadaas SP: 
Ameloblastin fusion protein enhances pulpal healing and dentin 
formation in porcine teeth. Calcif Tissue Int 2006, 2006(78):278-284. 

56. Tamburstuen MV, Reseland JE, Spahr A, Brookes SJ, Kvalheim G, Slaby I, 
Snead ML, Lyngstadaas SP: Ameloblastin expression and putative 
autoregulation in mesenchymal cells suggest a role in early bone 
formation and repair. Bone 201 1, 48:406-413. 

57. Turner P: Some observations on the dentition of the narwhal {Monodon 
monoceros). J Anat Physiol 1872, 7(Pt 1):75-79. 

58. Ishiyama M: Enamel structure in odontocete whales. Scanning Microsc 
1987, 1:1071-1079. 

59. Nweeia MT, Eichmiller FC, Hauschka PV, Tyler E, Mead JG, Potter CW, 
Angnatsiak DP, Richard PR, Orr JR, Black SR: Vestigial tooth anatomy and tusk 
nomenclature for Monodon monoceros. Anat Rec 201 2, 295:1 006-1 01 6. 

60. Zeichner-David M, Chen L-S, Hsu Z, Reyna J, Caton J, Bringas P: Amelogenin 
and ameloblastin show growth-factor like activity in periodontal 
ligament cells. Eur J Oral Sci 2006, 1 14(S up pi 1):244-253. 

61. Vymetal J, Slaby I, Spahr A, Vondrasek J, Lyngstadaas SP: Bioinformatic 
analysis and molecular modeling of human ameloblastin suggest a two- 
domain intrinsically unstructured calcium-binding protein. Eur J Oral Sci 
2008, 116:124-134. 

62. Deutsch D, Haze-Filderman A, Blumenfeld A, Dafni L, Leiser Y, Shay B, 
Gruenbaum-Cohen Y, Rosenfeld E, Fermon E, Zimmermann B, Haegewald S, 
Bernimoulin J-P, Taylor AL: Amelogenin, a major structural protein in 
mineralizing enamel, is also expressed in soft tissues: brain and cells of 
the hematopoietic system. Eur J Oral Sci 2006, 1 14(Suppl 1):183-189. 

63. Haze A, Taylor AL, Blumenfeld A, Rosenfeld E, Leiser Y, Dafni L, Shay B, 
Gruenbaum-Cohen Y, Fermon E, Haegewald S, Bernimoulin J-P, Deutsch D: 
Amelogenin expression in long bone and cartilage cells and in bone 
marrow progenitor cells. Anat Rec 2007, 290:455-460. 

64. Delgado S, Ishiyama Y, Sire J-Y: Validation of amelogenesis imperfecta 
inferred from amelogenin evolution. J Dent Res 2007, 86:326-330. 

65. Al-Hashimi N, Sire J-Y, Delgado S: Evolutionary analysis of mammalian 
enamelin, the largest enamel protein, supports a crucial role for the 32- 
kDa peptide and reveals selective adaptation in rodents and primates. 
J Mol Evol 2009, 69:635-656. 



doi:1 0.1 1 86/1 471 -21 48-1 3-20 

Cite this article as: Meredith et al.: Molecular decay of enamel matrix 
protein genes in turtles and other edentulous amniotes. BMC 

Evolutionary Biology 201 3 1 3:20. 



